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We show that the literature on pion exchange between charm and bottom mesons is inconsistent. 
We derive the formalism explicitly, expose differences between papers in the literature and clarify 
the implications. We show that the X(3872) can be a bound state but that results are very sensitive 
to a poorly constrained parameter. We confirm that bound states in the BB sector are possible. The 
circumstances whereby exotic combinations can bind with cc or bb quantum numbers are explored. 



I. INTRODUCTION 

The nature of the enigmatic charmonium meson A (3872), which appears at D°D*° threshold with J p = 1+, has 
been the subject of intense debate ever since its discovery. Historically Ericson and Karl[l| considered pion exchange 
in hadronic molecules. Tornqvist predicted 0, H[ that one pion exchange between charmed mesons gives an attractive 
force in the 1 = channel such that molecular or resonant DD* states might arise near threshold. Following the 
discovery of the A(3872) several papers suggested that it could be a 1 ++ state, driven by pion exchange HH and/or 
quark exchange^;] where flavour symmetry breaking was associated with its affinity for the neutral D°D*° threshold. 

More recently several papers have appeared assessing the potential attractive forces and asking whether a bound 
state is dynamically realisable. These notably include Suzuki0] who has argued that the one pion exchange forces 
are only able to make a feeble attraction at best, and most recently Liu et. al.[l, Q who have claimed that a bound 
state does not exist for reasonable values of parameters. 

All of these papers [1, [H, 0, 0, H, Q make different assumptions of detail, are not always self-consistent, and do not 
all agree on the mathematical expressions even where their assumptions are the same. Hence the purpose of the 
present paper is to attempt a unified treatment of this problem, enabling comparison between the various approaches 
to be made. In particular we shall make explicit the calculation of some critical signs, upon which attraction or 
repulsion can depend, and whose derivation is not described in the existing literature. As a result we shall find that 
expressions in the various papers are mutually incompatible. We shall then propose a consistent formulation, discuss 
its consequences and compare with the existing literature. 

Tornqvist @ initially assumed isospin symmetry and found 1 = attraction. Following the experimental discovery, 
Close and PageQ showed that the d — u mass difference can lead to substantial breaking of the flavour symmetry 
enabling attraction in the neutral D°D*° configuration. Tornqvist also studied isospin symmetry breaking[4j. Recently 
Liu et. al.[8, 9], using the empirical fact that the state is at neutral threshold, have focused solely on the neutral channel 
without discussion of isospin symmetry or its breaking. In addition, they assume that any phenomenon at the BB* 
threshold also involves only one charge channel. As we shall discuss here, the mechanism of flavour symmetry breaking 
can be critical in deciding which channels if any are attractive, and in particular whether attractive forces are strong 
enough to bind. We will also note that the BB* and DD* situations can be very different. As we shall argue, the 
symmetry breaking and dynamics in the DD* relative to BB* cases depend on the mass splittings between vector 
and pseudoscalar masses and whether they are larger or smaller than the ir mass. This has been noted clearly in the 
work of SuzukiQ and Liu et. al.[l,@ where the Fourier transform gives different potentials in position-space, but is 
not apparent in the original work of Tornqvist Q. Furthermore, we find differences in some critical signs relative to 
Tornqvist in Ref. [jj. These have potential implications for the attraction or repulsion in the DD* and BB* systems 
which differ from that reference. 

In the present paper we shall first derive the expression for the 7r-exchange potential along the lines of the original 
paper This will expose the origin of the signs that determine the overall attraction and repulsion and how the 
vector-pseudoscalar mass gap is critical. We shall concentrate on making contact with existing literature, showing 
where there are differences of assumption, sensitivity to inputs, and possible errors of calculation. Finally we shall 
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assess the implications. 

In Section III Al we make pedagogic comments about different conventions for charge conjugation eigenstates, in 
order to clarify discussions in the literature and to define our formalism. In Section III Bl we give a simple illustration 
of the spin expectation values for PV (VP) and VV. This exposes a relative sign between these that disagrees with 
Ref. 0. In Section EC] we give the overall spin and flavour factor. We calculate the effective potential in position 
space in Section [TTTl our result in equation (23) exposes the differences in the existing literature [3, 0, B @- 

In Section IIVI we move on to applications of the formalism: we discuss normalisation in Section IIV Al show the 
shape of the potentials in Section [IV Bl and apply the formalism to the DD* system in Section HV CI In Section [TV Dl 
we study the BB* and flavour exotic DD* and BB* systems. We finish with some general comments and conclusions 
in Section [Vl 



II. BASICS 

The effective one pion exchange potential between two light quarks can be split into a central term proportional to 
(<r, • 0j)(t, ■ tj) and a tensor term proportional to Sij(r)(Ti ■ Tj) with Sij(q) = 3(<r, • q)((Tj ■ q) — (<r, ■ (Tj). Ti are 
the isospin matrices acting on light quark i and cr^ are the Pauli spin matrices acting on light quark i. The matrix 
elements of these operators capture all the spin and flavour dependence. The remaining dependence on kinematics, 
normalisation of the potential and more detailed model assumptions is discussed in Sections IIIII and IIVI 

We obtain the effective potential between two hadrons by summing over all the interactions between light u and d 
quarks/ antiquarks. In heavy-light mesons such as D and B there is only one possible interaction to consider. 



A. Charge conjugation and conventions 

There has been some confusion in the literature as to the correct definition of the C = + state made of DD* ± c.c. 
In Ref. [3j the state with definite C parity is defined by (PV)± = [PV ± C(PV)]. As C 2 = 1 the above equation is 
self-consistent but does not specify the wavefunction until one chooses a convention whether CV(P) = ±V(P). Only 
the neutral state is an eigenstate of C and so one has CV = ±V = —V, which is consistent with CV = ±V = —V. 
The intermediate step CV(P) — ±V(P) is arbitrary (i.e. convention dependent). 

In Ref. 0, Q the C = + state is defined to be DD* + DD* , whereas in Refs. [j| B [l(| it is claimed that the 'correctly 
argued' form is DD* — DD* . As the overall sign of attraction versus repulsion depends on this sign, it is important 
to understand the origin of these alternate forms. As a particular example of our discussion above, the D and D* are 
not eigenstates of C, and so the eigenvalues for the eigenstates DD* ± DD* depend on what conventions are used to 
define the states. 

Hence first we define our flavour states. For a q and q of given flavour and spin at positions 1 and 2, eigenstates of 
C for P = CT+ and V = 1— are 

P:(q(l)q(2)+q(l)q(2))/V2 

V:(q(l)q(2)-q(l)q(2))/V2 (1) 
For a meson made of a heavy quark, Q, and light flavour q, we define the meson and 'anti-meson' to be 

P : (g(l)Q(2) + Q(l)?(2))/V2; P : (g(l)Q(2) + Q(l)g(2))/V2 

V : (g(l)Q(2) - Q(l)«(2))/V2; V : (g(l)Q(2) - Q(l)q(2))/V2 (2) 
With these definitions the C-eigenstates are 

C\PV±PV) =±\PV±PV). (3) 

As Tornqvist 0] used this convention, we shall do so in order to make most immediate comparison with his results. 
We agree however that there are possible advantages in using the opposite convention, for example, in making contact 
with approaches that use interpolating currents in QFT as argued in Refs. [§. l9l. Iic| . 
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B. Quark states and one pion exchange 

Our plan in this section is to calculate the sign of the spin expectation values of the central term in two cases: (i) 
PV or VP in J z = +1 and (ii) VV in J z = 2. We will show that the C = + combination of (PV + VP) has the same 
sign for the spin operator as does VV in J z — 2. This is opposite to the results in Tables 2 and 5 of Ref. 

As the 7T vertex connects D «-» D* (antiparticles understood here also) we will calculate 

(PV ± PV\H\PV ± PV) = 7](PV ± PV\H\VP ± VP) (4) 
so that the sign of the effective matrix element will be given by 

± V (P\h\V)(V\h\P) (5) 

The overall sign is therefore dependent on rj. 

From general arguments, we have ij = 1 for the J p — 1 + state: because L is even there is no phase change from 
the spatial wavefunction on swapping V(P) and P(V). There is no phase from the spin wavcfunction because spin 
1 coupling with spin to give spin 1 has the same sign as spin coupling with spin 1 to give spin 1. To show 
this explicitly we need to define the quark content of the mesons, including their spin orientation. For simplicity we 
consider the vector to be in state J z — +1, and for shorthand write [Qq] = Q^q^ (and similar for qQ etc), while the 



5 = state is denoted (Qq) = W \(Q* q l - Q [ q 
This implies that \PV) and \ VP) are as follows: 

\PV) = (Qq)[qQ] - (Qq)[Qq\ - (qQ)[Qq\ + (qQ)[qQ] (6) 

and 

\VP) = -lQq](qQ) - [Qq](Qq) + [qQ](Qq) + [qQ](qQ) (7) 

and thus i] = 1 in Equs.[4]and[5] 

It is sufficient to consider one ordering as long as we are consistent in our conventions. We choose A(Q 1 q 2 )B(Q s q4) — > 
A' '(Qiq' 2 )B' (Qsq 1 \) where q and q are the light u and d quarks and antiquarks, and Q and Q are the heavy quarks 
and antiquarks. In this convention the spin and isospin operators act on the light quarks/antiquarks i — 2 and j — 4. 

For simplicity of presentation just consider the terms where the heavy Q and Q are in positions 1 and 3 respectively; 
the analysis trivially applies to all other combinations with the same conclusion. The relevant terms are then 

\PV) = (-QU l Q ] q ] + Q l q ] Q ] q V ) (8) 

and 

\VP) = (-QU^Q^q 1 + Q\ ] Q l (?) (9) 

The 7r exchange leaves the spins of the Q(Q) unchanged and the V <-> P spin transition comes from the q(q). The 
non-zero transitions are then between the first term in Equ. [S]and the first term in Equ.O In each case the operator 
cr_(7 + gives +1 and hence the overall sign from the spin contributions to (PV\H\V P) = + 

Hence the sign is 

(PV ± PV\H\PV ± PV) = (PV ± PV\H\VP ± VP) =± = Sign(C) (10) 
The case of VV with all spins aligned, S = S z = 2, in the conventions above is 
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|FF) = [Qq] [Qq] - [Qq] [qQ] - [qQ] [Qq] + [qQ] [qQ] (11) 

Here again, focusing on the terms where the heavy Q and Q are in positions 1 and 3 respectively, and noting that the 
Q(Q) spins do not flip, it is immediately obvious that the sign is positive: 



(VV\H\VV}=+ (12) 

and hence the same as that for the (PV + PV), C = + channel 1 ++ . 

The explicit inclusion of flavour (isospin) for the qq = dd±uu (and appropriate charge conjugated form) introduces 
further signs, causing I — and 7=1 channels to have opposite behaviours. However these factors are common to all 
of the above and do not change the general conclusion that for a given isospin the C — + combination of (PV + VP) 
has the same sign as does VV in J z = 2. 

It is possible to deduce the overall sign as follows. 

Start from TV TV with 5 = 1 and 7 = 0: the deuteron. This involves it exchange between q and q; there is no q — q 
interaction here. Now consider the case of VV with 5 = 2 and 7 = 0, which is like the deuteron in that again there 
is no q — q interaction. The (<Xj • (Tj) has the same sign in both cases. The tt exchange (G-parity) gives opposite sign 
between q — q and q — q, and hence VV with 5 = 2 and 1 = has opposite overall sign. So far everything agrees with 
the calculations in Ref. [H. It is only now, where Equs. fTOl and fT2l imply that the C = + combination of (PV + VP) 
has the same sign as does VV in 5 5 2 , and hence opposite to the deuteron, in contrast to Ref. [3]. 

We shall show that when the spatial matrix elements are calculated, we agree with the formulation in Refs. @, H, 
and disagree with Q . This ironically introduces a further relative sign in the spatial contribution of Ref. [|[ in the 
charm sector, which, as we shall show, will eventually cause the DD* + DD* state to be mildly attractive, in part as 
a result of two sign errors mutually cancelling. However, the spatial sign-flip does not occur in the heavy quark limit, 
and hence some care is required in comparing BB* + BB* to DD* + DD*. 



C. Overall Sign 



The overall sign is determined, inter alia, by the expectation values of (<Ti ■ tTj)(Ti ■ Tj) and Sij(q)(Ti ■ Tj). 

The spin matrix element of the central term, (cTi - (Tj), can be calculated explicitly as above or using general angular 
momentum theory. General expressions are given in Appendix [XJ It is straightforward to show that the spin matrix 
elements for PP -> PP, PV -> PV and VP -> VP vanish. For PV -> VP and VP -> PV we get +1 and for 
VV VV we get ^(St(St + 1) — 4). These results agree with the calculation given above (Section [II Bp . 

The spin matrix element of the tensor term, Sij (q) = 3(<Tj • q)(<Tj ■ q) — (&i • &j), can be calculated in a similar way, 
see Appendix [A] 

The isospin factor is trivial to evaluate: 

(T i -T j ) = ±((T i +T j f-T?-T?). (13) 

For two isospin half mesons interacting, this is 21(1 + 1) — 3, i.e. —3 in total isospin I = or +1 in I = 1. 

The interaction potential between a quark and an antiquark is opposite to that between two quarks (or two 
antiquarks) because of the G-parity of the pion. 

Adopting the above results, for total spin S, isospin states with isospin / and charge conjugation parity C, the 
flavour and spin factors for the central term are: 

• VV: (5(5 + 1) - 4)(7(7 + 1) - 3/2) 

• VV: -(5(5 + 1) - 4)(7(7 + 1) - 3/2) 

• PV: (27(7 + 1) - 3) 

• PV: -(7(27(7 + 1) - 3) 
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The VV and VV expressions agree with those of TornqvistQ. However, there are different overall minus signs in 
the PV and PV expressions compared to those of Tornqvist. 

For reference we note that the matrix elements for DD* with J p = 1 + , charge conjugation parity C and isospin I, 
in the basis L = 0, L = 2 are 



C {21(7+1) -3} 



1 
1 



V c 







-V2 



-V2 1 



(14) 



where Vc and Vp distinguish the central and tensor terms. 

When isospin symmetry is broken, we replace the (t± ■ Ti) factor by a 2 by 2 matrix 

-1 -2 

in the basis of charged/neutral states 00, + — . In the isospin limit, it is easy to see that this has eigenvalues —3 and 
1 with the respective isospin eigenvectors. Close and PageQ and Tornqvist Q discuss isospin symmetry breaking. 

If one only considers the particular interaction D°D*° — > D*°D° and not every charge combination, the interaction 
strength is reduced by a factor of 1/3 compared to the isospin symmetry limit with 7 = 0. Liu et. al.[8j only consider 
this particular charge interaction. 

We can not really say whether the potentials arc attractive or repulsive until we know something about the kinematic 
dependence which we discuss in the next section. 



III. EFFECTIVE POTENTIAL IN POSITION SPACE 



If the interaction between n and a light quark is taken to be 

£ = ^-q(x)ri5Tg(x) ■ d M <t>(x), (15) 

JTT 

the effective potential between two light quarks due to one pion exchange in the static limit is 

v( . _ 9 2 fa • q)(gj- • q) , . rifi N 
f4 t - m i 

Here m w is the ir mass, / w the n decay constant (defined by V2 (0 |yLj(0)|7r° (q)) = if^q^ or (0|A /J (0)|7r ± ((7)) = if^q^ 
as in the PDG Review 2006[12]), g a dimensionless coupling constant, and q is the four-momentum transfer. For 
elastic scattering this reduces to the form given in Ericson and Weise[ll|: 

v(n\- ff 2 (ovgXgj-g) / n7 N 
v{q) - -72 — j — — ; — 2 — i Ti ' T i> y 1 ') 

fl \q\ 2 + ml 

V(r) = ^(n- TiHai ■ V)( CTj • V)^^ (18) 

where q is the three-momentum transfer. 

Following Tornqvist we define /j 2 = m\ — (my — nip) 2 , my is the vector meson mass and mp the pseudoscalar 
meson mass. Close to the static limit (both initial and final 3-momenta are zero) we have g§ « [my — mp) 2 and so 
q 2 — m 2 sa — \q\ 2 — /i 2 . In the limit of elastic scattering /1 = m n . 

We write the overall scale in terms of Vb : 

« s ^ (19) 

and so the expression for the effective potential is 

v( n \- mkVq (v* ■ q)(<rj ■ q) r _ t) (9n) 
y(9) -"^| \ q] 2 + fl2 in -Til (20) 
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This is the same as Tornqvist's except for an overall minus sign - this missing sign in his momentum space form turns 
out not to be important because his expression for the potential in position space does have the correct sign. 

Note that there are other definitions of (j, in the literature. SuzukiQ and Liu et. al.Q use 2(mv 



nip 



—fJr, where the correspondence, with a minus sign, is exact in the static limit. As mentioned by Suzuki, the form of 
the potential in position space depends on the sign of /z 2 . 

Following Ericson and Weise[ll[, it is useful to highlight the different spin dependences in the potential. The central 
potential is proportional to (<Tj • crj) (hyperfine interaction-like) and the tensor term is proportional to Sij(q): 



V(q) = 



AttV 



-1 



(a-i -crj)- 



(21) 



with Sij(q) = 3(<Tj • q)(cTj ■ q) — (er,; • crj). The literature 0, d, 0, H, 0> H, HH nas focused on the potential in position 
space but disagree on its form. Hence we take the Fourier transform and compare with the literature. 

For ji 2 > 0: 



V(r) = V Q 



Ait 
mi 



fi 2 e-^ r 
m 2 m^r 



[i 2 e'^ 
m 2 m^r 



3 

1 + — 

[IT 



fj, 2 r 2 



Sij(r) 



and for /i 2 = — jl 2 < the real part is 
V(r) = V 



47T 



jl 2 cos jlr 
m 2 m^r 



(cr, • crj) - 



M 2 



m,±. m^r 



cos [ir 



3 sin fir 3 cos fir 



fir 



\l 2 r 2 



S l3 (r) 



(22) 



(t 4 • Tj ) (23) 



where we have kept the terms in the same order. 

The expression with /i 2 > reproduces that of Tornqvist 3j, except that he ignores the 5(r) term. He allows for a 
varying /j, 2 but his quoted potential appears implicitly to assume that /j, 2 is always positive. In cases where /j, 2 < 0, our 
potential agrees with that of Liu et. al.@| (except that they have ignored the tensor term) and with that of Suzuki Q 
apart from some relative minus signs (where we confirm the form of Liu et. al.). 



Regularization 



The potential is singular at small distances and so has to be regularised. Following Tornqvist and Liu et. al., we do 
this by introducing a form factor at each n vertex which leads to an extra factor of ^ \i m q % ^ multiplying V(q) (Equ. 



[21]) . This gives the n an effective RMS radius of y/To/A3\. The value of A has to be determined phenomenologically. 
Tornqvist mentions that in nucleon-nucleon interactions values between 0.8 and 1.5 GeV have been used depending 
on the model and application, but that larger values (A > 1.4 GeV) are required for nucleon-nucleon phase shifts. 
He says that for heavy mesons which have a smaller size than nucleons, one would expect a smaller effective radius 
of the 7r source corresponding to a larger A. A larger A gives a stronger potential at short distances; we shall find in 
Section HVl that the results depend strongly on A. 



The central terms (i.e. all but the tensor Sij(r)) for /s 2 > are 



A(A 2 -m 2 ) t 
2mi 



-Xr 



-fir 



-Xr\ 



(24) 



with X 2 = A 2 



Tornqvist omits the l S function' piece (more precisely the piece that gave the delta function 



when no regularisation was used), "which from the phenomenological point of view will be included in the short range 
potential and regularisation scheme" |3J. This results in all his terms having a common /i 2 factor. If for comparison 
we also ignore the l 5 function' piece we get 



Vr, 



l/i 2 (A 2 
2 



Xmi 



m l) e -Xr 



-Xr 



(a-i ■Oj-XTi - Tj) 



Note that this is not quite the same as the expression in Ref. namely: 



(25) 



M 2 (A 2 - M 2 ) c - Ar 
Ami 



,— Ar 



((Ti ■ CT 3 ) {n -Tj). 



(26) 
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The most significant difference (e.g. see Figures [T] and [2]) is due to the missing factor of 1/2 in the first term. The 
other differences are a m 2 — > fi 2 in one place and the approximation X f» A. 

We now comment on the case /x 2 = — ji 2 < 0. The real part of the full central term (including the '5 function' piece) 

is 



X(A 2 -ml) Xr p, 2 , x . 
3 - e- Xr - (cos H - e- Xr ) 



(<7; '(Tj){Ti -Tj). 



This expression agrees with Equ. 19 of Liu et. al. If we ignore the l S function' piece we get 

/i 2 (A 2 -m 2 ) Xr £ 2 



2Xmi 



(cos(/tr) - e~ Xr ) 



fa -cr^in - Tj), 



The tensor terms for fi 2 > are 
V(r) = V 



,,2 p-ttr / o 



TOt Vfl-rrT 



(27) 



(28) 



(29) 



X 2 e- Xr ( . 3 



m 2 m^r 



Xr X 2 ? 



(A 2 - ml) e~ Xr 
2m 2 m^r 



(1 + Xr) 



S ij(r){l~t ■ Tj 



This agrees with the Equ. 26 of Tornqvist apart from one place in the last term where m 2 — » /x 2 . For /i 2 = — /x 2 < 
the real part of the tensor terms are 



V(r) = V 



m 2 m^r 

2 a -Xr 



cos fir 



3 sin fir 3 cos fir 



fir 



fi 2 r 2 



(30) 



X z e 



mtL m^r 



Xr X 2 r 2 



-Xr 



(A 2 - m 2 ) e 
2m 2 m x r 



(1 + Xr) 



SiJf^Ti-Tj). 



We collate some useful relationships and Fourier transforms in Appendix [B] from which the derivation of the above 
expressions may be checked. 

We shall use the forms Equs. [2^1 and 1 2 71 for the central term, and Equs.l2"9"landl3"0"lfor the tensor term in our analysis. 
We shall illustrate the implications for binding by comparing with other expressions in the previous literature. 



IV. APPLICATIONS 



Normalisation 



We have defined Vq in the same way as Tornqvist [3j and so we can write it in terms of the 7r-nucleon coupling 
constant f^N'- 



Air 



25 « 2 < = 25 3^ ^ 08 



9 / 2 4tt 9 m, 

Tornqvist finds Vq ~ 1.3 MeV, which is consistent with Ericson and Weise[ll[ and Ericson and Karl[l|. 
Vo can also be related to the D* — > D°ir + transition: 

T(D* - D\+) = -fi- P l = 2V ^. 

6nf£ m% 



(31) 



(32) 



The PDG Review 2006[T3| gives T(D* D°tt+) = (65 ± 15)keV and so V = (1.5 ± 0.3) MeV which is consistent 
with the above. 

Suzuki[3] appears to have an effective Vq ~ 0.73 MeV and Liu et. al.Q have Vq ~ 0.68 MeV. Note that there are 
various definitions of the n decay constant used in the literature. We take the definition of the PDG Review 2006 (l2j 

where « 130 MeV. Compare this with the definition used in SwansonQ where — > At and = 92 MeV. 
We will take Vq = 1-3 MeV throughout this work. 
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B. Shape of The Potentials 

In Figures [1] and [5] we plot the central part of the potential PV in I = 0, L = 0, J PC = 1 ++ with A = 1 GeV, 
Vo = 1.3 MeV, and fj, 2 = m 2 or /i 2 = —2000 MeV 2 respectively. As well as our expression for the full central 
potential (Equs. |2~H and |2~T|) , for comparison we plot the potential without the '5 function' piece (Equs. [2"5l and |2"5|) 
and Tornqvist's expression for the potential (Equ. I26p . Note that in all three cases we use our expression for the spin 
and flavour factors as discussed in Section Hi CI 

Figure Q] shows that for positive fj, 2 the '5 function' term is attractive, dominates at short distances and is opposite 
in sign to the other central terms. The effect of the different factor of 2 in Tornqvist's potential compared to ours 
when we ignore the '<5 function' term can also be seen, in particular, the sign change as r — > 0. The three potentials 
differ at short distances but become indistinguishable for r > 2 fm. 

Figure [2] shows that for negative /z 2 the '5 function' term remains attractive, dominates at short distances and is 
the same sign as the other central terms. The effect of the different factor of 2 in Tornqvist's potential can again 
be seen, as can the fact that the three potentials differ at short distances but are the same at long distances. As is 
apparent from comparing Figures [T] and [2l the '6 function' term does not change sign when /i 2 — > — fi 2 whereas the 
other central terms do. Thus the phenomenological conclusions may depend significantly on how the 5-function term 
is treated. 

In Figure[3]we plot our expression for the central potential with the above parameters, fi 2 — —2000 MeV 2 and show 
the effect of varying A. It can be seen that increasing A increases the strength of the potential at short distances and 
that the potential is quite sensitive to A. 



V(r) / MeV 




FIG. 1: Central potential for PV in I = 0, L = 0, J PC = 1 ++ with A = 1 GeV, V = 1.3 MeV, and /i 2 = ml. (1) is our 
potential (Equ.[2U), (2) is our potential without the '8 function' piece (Equ. [231) . to be compared with Tornqvist's analogous 
(i.e. no Vfunction) potential (Equ. [26} shown as curve (3). 



C. Application to DD* 

In this section we discuss application of the above formalism to the DD* system with J PC = 1 ++ . The mesons 
can be in two different partial waves (L = or 2) and mixing between these may be important. 



9 



V(r) / MeV 




FIG. 2: Central potential for PV in I = 0, L = 0, J PC = 1 ++ with A = 1 GeV, Vo = 1.3 MeV, and fi 2 = -2000 MeV 2 . (1) is 
our potential (Equ. 127). (2) is our potential without the '6 function' piece (Equ. I28|) . (3) is Tornqvist's potential modified for 
negative fj 2 . 



We solve the potential by discretising the time-independent Schrodinger equation with the one pion exchange 
potential and then diagonalising the resulting matrix. We vary the number of points and the maximum radius to 
check for discretisation and finite volume effects. 

To get some handle on the size of A, we first calculate and solve the effective potential for the deuteron. We 
assume isospin symmetry and include both the central and tensor terms. For comparison, we consider three different 
expressions for the central potential: (1) our full expression (Equs. [Ml and [77]) . (2) Our expression without the l S 
function' piece fEqu. I2"5l and 125)) . and (3) Tornqvist's expression (Equ. 123)) . We vary A to fix the binding energy to 
w 2.2 MeV. 

For potentials (1), (2) and (3) we find the values required are A = 960 MeV, 750 MeV and 760 MeV respectively. 
More details of the calculations and results are given in Appendix [C) A gives an effective radius of the pion ~ 1/ A 
which is expected[3j to depend on the interacting hadron and so be smaller for the DD* system compared to the 
deuteron. This argument implies that the A we have found from solving the deuteron is a lower bound on the value 
of A that should be used for the DD* system. The value of A required for the BB* system should be larger again. 

Moving to the DD* system, we verify the results of Ref. 8] within the assumptions made there (noting the 
possible difference in normalisation discussed above in Section IIV A[) . Specifically, we initially considered only the 
D°D*° + D°D*° charge combination, used our full expression for the central term (1) and ignored the tensor term. In 
this case the flavour factor is —1 and one pion exchange is slightly inelastic with a negative /i 2 {p 2 = —2000 MeV 2 ). 
For a reasonable range of A (we checked up to 2.5 GeV) we find no bound state, supporting the results of Ref. Q 
within their assumptions. However, for a larger A of 3 GeV, we do find a bound state with energy sa 3868 MeV. 

However, this shows how sensitive the results are to such assumptions. For example, if we now assume isospin 
symmetry with 1 = (flavour factor is —3) and don't change anything else, we do find a bound state. For A in the 
range 960 MeV to 1200 MeV we find E = 3871 MeV - 3863MeV. This highlights the importance of including all the 
possible charge combinations. 

We then performed a more complete calculation including the tensor term. We also allow for isospin symmetry 
breaking through different charged/neutral meson masses and [i 2 s. The ji 2 s are all negative except for D*° — > D + 
which has a positive /j, 2 . Because we can only use one /i 2 for each charge mode, there is ambiguity as to whether we 
should calculate the charged-neutral /i 2 from (A) D* + — ► D° or (B) D*° — > D + . We will compare the results obtained 
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FIG. 3: Central potential for PV in I = 0, L = 0, J PC = 1 ++ with Vo = 1.3 MeV, fi 2 = -2000 MeV 2 , and with our expression 
for the potential (Equ. I27|l for various values of A. 



with both choices (A) and (B). 

The detailed results are given in Table Hill of Appendix [Dl With our full expression for the potential (1), we find a 
bound state with E = 3870 MeV - 3806 MeV for A » 960 MeV - 1400 MeV. For A = 800 MeV we find no bound 
state. This shows that the binding energy is very sensitive to A. The choice of fj, 2 (A) or (B) modifies the binding 
energy by at most 1 MeV. In particular, Table IrU shows that isospin symmetry breaking is important, especially when 
the state is close to threshold. This is not particularly surprising: the charged-neutral mass difference gets relatively 
more important when the state is only just bound. The state is mostly L = with a small L — 2 component. 

Using potentials (2) and (3), we require a larger A for a bound state (A > 1750 MeV) because these potentials do not 
have the strong attraction of the l S function' term. We find E = 3871 MeV- 3863 MeV for A w 1750 MeV- 2000 MeV. 
The differences between potentials (2) and (3) do not significantly change the binding energy or other properties, as 
shown in Table [TTTJ 

In summary, we have found a bound state for reasonable values of A using all three potentials. The binding energy 
is sensitive to the potential and the value of A used. Clearly the exact values of the binding energies calculated 
are not significant, because of their sensitivity to parameters and the potential. Nonetheless, the general conclusion 
that within one pion exchange potentials a bound state can be formed appears to be robust. Conversely, there is 
no theoretical reason to expect that it must be formed, nor that it cannot occur. Basically, the phenomena of the 
V(3872), within the assumption of one 7r-exchange, would constrain the parameter A more than present knowledge 
of this parameter can constrain the DD* dynamics. 



D. Further Applications 



We can also apply a similar analysis to the analogous BB* system. Furthermore isospin breaking is expected to be 
less important here because of the relatively smaller charged/neutral differences. The larger B meson masses mean 
that the kinetic energy is less of a hindrance to binding. 

Results are given in Table IIVI of Appendix [D] which show that a bound state is found for reasonable values of A, 
the binding energy is sensitive to A and isospin breaking is less important than in the DD* system. We note that 
there are larger L = 2 components here compared to the DD* system, especially when using potentials (2) and (3). 
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Because /i 2 is positive here, the central terms other than the '6 function' term are repulsive and so the tensor term is 
important for forming a bound state. 

The same analysis again can be applied to flavour exotics, that is states such as DD* and BB* with two charm or 
bottom quarks. The one pion exchange interaction in such states has opposite overall sign to that in the DD* and 
BB* systems (Section IHCj) . However, mixing with the tensor terms complicates the situation and we must do more 
than just consider the sign of the central term. 

For the DD* system, results are given in Table fVl of Appendix [Dl With potential (1) we find a bound state with 
energy 3870 MeV for A = 2000 MeV. There is no bound state for A = 1700MeV. The importance of the tensor 
term is apparent from the relatively large L — 2 component. The +0 and 0+ components have the opposite sign and 
so the state would be isoscalar in the isospin symmetry limit. Using potential (2) or (3) we find a bound state with 
energy w 3870 MeV for A = 1500 MeV, a bound state with energy w 3810 MeV for A = 2000 MeV, and no bound 
state for A = 1000 MeV. Again, there are relatively large L = 2 components. 

For the BB* system, results are given in Table IVTl of Appendix iDl With potential (1), we find a bound state with 
energy w 10600 MeV with A = 1000 MeV, a bound state with energy w 10560 MeV with A = 1500 MeV, and no 
bound state for A = 700 MeV. There are again relatively large L = 2 components and the +0 and 0+ components 
have opposite sign. Using potential (2) or (3) we find a bound state with energy w 10600 MeV for A = 700 MeV, a 
bound state with energy w 10590 MeV for A = 1000 MeV, and no bound state for A = 500 MeV. Again there are 
relatively large L = 2 components. 

In summary, we find great sensitivity to parameters, but qualitatively confirm that BB* and exotic states can bind 
in one pion exchange. As can be seen from Table U using our expression for the potential (1), we find DD*, BB* 
and BB* can bind with A ~ 1000 MeV. However, DD* requires a larger A ~ 2000 MeV. Hence, within our potential 
and assumptions, the parameters that allow V(3872) to emerge as a bound state preclude binding the exotic DD* 
channel. However, binding in both BB* and exotic BB* are possible. 



System 


Potential 


Approximate A/MeV 


DD* 


1 


960 




2 and 3 


1750 


BB* 


1 


< 800 




2 and 3 


1000 


DD* 


1 


2000 




2 and 3 


1500 


BB* 


1 


1000 




2 and 3 


700 



TABLE I: Approximate minimum values of A required to bind 



V. CONCLUSIONS 

In summary, we agree with the qualitative results of Tornqvist, but as a result of various differences cancelling out. 
SwansonQ found that quark exchange alone did not bind within the one gluon exchange contact approximation. This 
led him to include one pion exchange based upon the work of Ref . • 

We have quantified the arguments of Ref. 0] that a small fi 2 leads to a small binding energy. In doing so, we have 
included the tensor term, and also flavour factors, neither of which were discussed in that reference. Our calculations 
show that the sensitivity to A is the overriding factor. 

Liu et. al.[8j considered only D°D *; charged modes were ignored, as was the tensor term. Within their assumptions 
we confirm their results, though there is the question of overall normalisation. Our work highlights the importance 
of taking into account all charged modes. They have also considered a exchange, and argue that this makes it harder 
for the DD* to bind. 

We have discussed flavour exotic states and find that the overall sign of the central term alone does not determine 
whether or not a bound state is formed; mixing due to the tensor term is also important. It can be dangerous to 
ignore the tensor term and its contribution to different processes can be important: this is especially true for the case 
of flavour exotics. 
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APPENDIX A: SPIN MATRIX ELEMENTS 



We choose A(Qiq2)B(Q 3 q~4) — + A' (Qiq 2 )B' (Q 3 q' 4 ) where q and q are the light u and d quarks and antiquarks, and 
Q and Q are the heavy quarks and antiquarks. In this convention the spin operators act on the light quarks/antiquarks 
i = 2 and j = 4. 

In general, for initial state mesons with spin Sa and Sb coupled to total spin St going to final state mesons S' A 
and S' B coupled to total spin S T , the matrix element of the central term, (<T2 ■ C4) is given by: 

{S' A S' B {S' T M')\{(T2 ■ a 4 )\S A S B (S T M)) = (Al) 

1/2 1/2 S A 1 f 1/2 1/2 S 13 
2(S 24 (S 2 4 + l)-3/2)^s A s B s' A s' B s 13 s 13 s 2i s 2i < 1/2 1/2 S B \{ 1/2 1/2 S 2i 



$s' t ,s t 8m',m / J 



S13S24 



5*13 5*24 S T J [ S' A S' B S' T 



where {} are Wigner 9j symbols and Uabc... = y/(2A + 1)(2B + 1)(2C + 1)... . 

From this expression, it is straightforward to show that the spin matrix elements for PP — > PP, PV — > PV and 
VP -► VP vanish. For PV -> VP and VP -> PV we get +1 and for VV VV we get ±(S T (S T + 1) - 4). These 
results agree with the direct calculation given above. 

The matrix element of the tensor term, S*24(g) = 3(<T2 • <z)(f4 • q) — (&2 • C4), can be calculated in a similar wav(l3j. 
With the spins defined above and with the total spin Jab, J'abi relative orbital angular momentum L, L 1 and total 
angular momentum J , J 1 = J in the initial and final states respectively, the matrix element is: 

(A'B'\S 2i (q)\AB) = (A2) 

JS13S24 

1/2 1/2 S A } [ 1/2 1/2 S' A ) r ~ 

1/2 1/2 S B 1/2 1/2 S' B \ \ g 13 J J AB >^ "V v )\ c ^ \ „: 2. (,) 

5*13 S'24 Jab J [ S13 S24 J'ab J 

where {} are Wigner 6j and 9j symbols and (L\ 0\L, 0; 2, 0) is a Clebsch Gordan coefficient 



SaSbS'a S' f ,Si3Si3S24S24S24.Jab J'a rJJL 



V 


S24 


f ) 


\ L ' 


L 


5} 


Sl3 


J 


Jab J 


y S24 


5*24 
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APPENDIX B: USEFUL EXPRESSIONS FOR CALCULATING THE POTENTIAL IN POSITION SPACE 



We collate useful Fourier transforms used in Table HT1 where 

1 



V(r) = 



(2k) 



V{q)e l i r d 3 q. 



(Bl) 



A useful decomposition is 

(o-i ■ g)(cr 2 ■ q) 

|q| 2 ±M 2 

with S 12 (q) = 3(<Ti ■ q)(er 2 ■ q) - (<Xi • er 2 
The form factor is \F(\q 



\q\ 2 +X'< 



(<ri . cr 2 j =F (ci • cr 2 j 



2M2 _ / A^-mJ ^ 



S 12 (q) 



(B2) 



with X 2 = A A + [i 1 — m 2 = A" — /Er — m~ and p, 



V(q) 



V{r) 



S(r) 



4?r r 
1 e^ r 
47T r 



M 2 +M 2 



2tt^ Jo 



Sll(q)j$^ 



-S 12 (f)£^(l + A + _ 



_3 



|f(k| 2 )| 5 



(A 



°° (|g|")l^o(|qk) 



2-k 2 Jo 



d\q\ 



i 

1 

477 



£X) -Xr 
(A 2 — m 2 ) „-Xr" 



-5ia(# 



^ l 

4tt 



— + 3 

/Jr (fir) 2 

Zi 3 
p.r 



-^la^J^Jo M 2 + M 2 d l9l 



(Ar)' 



1 



1 + 



.Y7- 



3 



1 + Xr + (Ir)! 



-(A 2 -m^)V(l+Xr) 
-(A 2 -m^(l+Xr) 



TABLE II: Summary of Fourier transforms used 



APPENDIX C: DEUTERON POTENTIAL 

The deuteron is a combination of two nucleons in an isosinglet state with J — 1, S = 1, L = or 2. Because we 
have to sum over interactions between all the light quarks, relative to the heavy-light mesons the spin factors are 

changed to — 2£ [ ^ ^ | for the central term and I ^- ^ I for the tensor termS. The mixing between the 

3 \^ o i y 3 \V8 -2 J M 

L = and L = 2 states is important in binding the deuteron. 

If we follow Tornqvist and assume that pion exchange is the only binding mechanism, we can require the binding 
energy to be ~ 2.22 MeV and so fix the scale A. We take Vq = 1.3 MeV throughout. 

We use three different expressions for the potential: 

1. Our full expression given in Equ. 1241 

2. Our expression without the l 5 function' piece, Equ. [25] 

3. Tornqvist's expression^ 
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Potential 


Set 


A/MeV 


E/MeV 


00, L = 


+-,L = 


00, L = 2 


+-,L = 2 


RMS Radius/MeV" 1 


(1) 


(A) and (B) 


800 


(not bound) 












(1) 


(A) 


962 


3870 


0.76 


0.22 


0.01 


0.01 


0.02 


(1) 


(B) 


962 


3870 


0.79 


0.19 


0.01 


0.01 


0.02 


(1) 


(A) 


1000 


3868 


0.69 


0.28 


0.02 


0.02 


0.01 


(1) 


(B) 


1000 


3869 


0.70 


0.27 


0.02 


0.02 


0.01 


(1) 


(A) 


1400 


3806 


0.50 


0.46 


0.02 


0.02 


0.003 


(1) 


(B) 


1400 


3807 


0.50 


0.46 


0.02 


0.02 


0.003 


(2) 


(A) and (B) 


752 


(not bound) 












(2) 


(A) and (B) 


1500 


(not bound) 












(2) 


(A) 


1750 


3871 


0.86 


0.11 


0.02 


0.02 


0.03 


(2) 


(B) 


1750 


3871 


0.90 


0.08 


0.01 


0.02 


0.03 


(2) 


(A) 


1800 


3870 


0.78 


0.17 


0.02 


0.03 


0.02 


(2) 


(B) 


1800 


3870 


0.81 


0.14 


0.02 


0.02 


0.02 


(2) 


(A) 


2000 


3863 


0.59 


0.33 


0.04 


0.04 


0.007 


(2) 


(B) 


2000 


3863 


0.60 


0.32 


0.04 


0.04 


0.008 


(3) 


(A) and (B) 


f 760 


(not bound) 












(3) 


(A) and (B) 


1500 


(not bound) 












(3) 


(A) 


1750 


3871 


0.86 


0.11 


0.02 


0.02 


0.03 


(3) 


(B) 


1750 


3871 


0.90 


0.08 


0.01 


0.01 


0.03 


(3) 


(A) 


2000 


3863 


0.59 


0.33 


0.04 


0.04 


0.007 


(3) 


(B) 


2000 


3863 


0.60 


0.32 


0.04 


0.04 


0.008 



TABLE III: DD* Results. The amplitudes squared of the different components are given along with the RMS radius of the 
00, L — component. The H — and 00 L — components have the same sign (i.e. would be an isoscalar state in the isospin 
symmetry limit). 



We solve the potential by discretising the time-independent Schrodinger equation and then diagonalising the re- 
sulting matrix. We take N = 500 or TV = 1000 points for each L and use maximum radii i?o = 0.10 MeV -1 and 
0.15 MeV" 1 to check for finite volume effects. We set fi = = 135 MeV and the reduced mass \i = l/(l/m n + l/m p ) 
with m n = 939.57 MeV and m p = 938.27 MeV. 

Using potential (1) we obtain a binding energy E = 2.23 MeV with A = 962 MeV and we find a relative amplitude 
squared of 0.93 in L = and 0.07 in L = 2 with a RMS radius in L = of 0.02 MeV" 1 . 

Using potential (2) we obtain a binding energy E = 2.20 MeV with A = 752 MeV and we find a relative amplitude 
squared of 0.94 in L = and 0.06 in L = 2 with a RMS radius in L = of 0.02 MeV" 1 . 

Using potential (3) we obtain a binding energy E = 2.22 MeV with A = 760 MeV and we find a relative amplitude 
squared of 0.94 in L = and 0.06 in L = 2 with a RMS radius in L = of 0.02 MeV" 1 . 

Although the results are sensitive to A, once this is fixed the other properties do not depend strongly on the details 
of the potential. 



APPENDIX D: TABLES OF RESULTS 
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Potential 


A/MeV 


£/MeV 


00, L = 


+-,L = 


00, L = 2 


+ -,L = 2 


RMS Radius/MeV -1 


(1) 


800 


10580 


0.47 


0.48 


0.03 


0.03 


0.004 


(1) 


962 


10540 


0.47 


0.48 


0.03 


0.03 


0.002 


(1) 


1000 


10530 


0.47 


0.48 


0.03 


0.03 


0.002 


(1) 


1400 


10270 


0.47 


0.47 


0.03 


0.03 


0.001 


(2) 


752 


(not bound) 












(2) 


1000 


10600 


0.40 


0.50 


0.06 


0.05 


0.010 


(2) 


1400 


10580 


0.41 


0.42 


0.08 


0.08 


0.004 


(2) 


1750 


10520 


0.41 


0.41 


0.09 


0.09 


0.003 


(3) 


760 


(not bound) 












(3) 


1000 


10600 


0.40 


0.49 


0.06 


0.05 


0.010 


(3) 


1400 


10580 


0.41 


0.42 


0.08 


0.08 


0.004 


(3) 


1750 


10520 


0.41 


0.41 


0.09 


0.09 


0.003 



TABLE IV: BB* Results. The amplitudes squared of the different components are given along with the RMS radius of the 
00, L = component. The H — and 00 L — components have the same sign (i.e. would be an isoscalar state in the isospin 
symmetry limit). 



Potential 


A/MeV 


_E/MeV 


0+,L = 


+0,L = 


0+,L = 2 


+0,L = 2 


(1) 


1700 


(not bound) 










(1) 


2000 


3869 


0.46 


0.36 


0.09 


0.09 


(2) 


1000 


(not bound) 










(2) 


1500 


3871 


0.49 


0.39 


0.06 


0.06 


(2) 


2000 


3814 


0.39 


0.38 


0.12 


0.12 


(3) 


1000 


(not bound) 










(3) 


1500 


3871 


0.49 


0.39 


0.06 


0.06 


(3) 


2000 


3814 


0.39 


0.38 


0.12 


0.12 



TABLE V: Exotic DD* Results. The amplitudes squared of the different components are given. The 0+ and +0 L = 
components have the opposite sign (i.e. would be an isoscalar state in the isospin symmetry limit). 



Potential 


A/MeV 


£/MeV 


0+,L = 


+0,L = 


0+,L = 2 


+0,L = 2 


(1) 


700 


(not bound) 










(1) 


1000 


10600 


0.33 


0.38 


0.15 


0.15 


(1) 


1500 


10560 


0.27 


0.28 


0.23 


0.23 


(2) 


500 


(not bound) 










(2) 


700 


10600 


0.38 


0.46 


0.08 


0.08 


(2) 


1000 


10590 


0.35 


0.36 


0.15 


0.15 


(3) 


500 


(not bound) 










(3) 


700 


10600 


0.38 


0.47 


0.08 


0.08 


(3) 


1000 


10590 


0.35 


0.36 


0.15 


0.15 



TABLE VI: Exotic BB* Results. The amplitudes squared of the different components are given. The 0+ and +0 L — 
components have the opposite sign (i.e. would be an isoscalar state in the isospin symmetry limit). 



